One-Electron Physics of the Actinides 
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We present a detailed analysis of the one-electron physics of the actinides. Various LMTO basis 
sets are analyzed in order to determine a robust bare Hamiltonian for the actinides. The hybridiza- 
tion between /- an spd- states is compared with the f — f hopping in order to understand the 
Anderson-like and Hubbard-like contributions to itineracy in the actinides. We show that both 
contributions decrease strongly as one move from the light actinides to the heavy actinides, while 
the Anderson-like contribution dominates in all cases. A real-space analysis of the band structure 
shows that nearest-neighbor hopping dominates the physics in these materials. Finally, we discuss 
the implications of our results to the delocalization transition as function of atomic number across 
the actinide series. 

PACS numbers: 71.30.-(-h, 71.27.-|-a, 71.15.Mb, 71.20.-b 



I. INTRODUCTION 
A. Background of the actinides 

It is well accepted^ that the actinides are divided into 
two groups based on the behavior of the /-electrons. The 
hghter actinides (Th to Pu) have smaller atomic volumes, 
low-symmetry crystal structures and itinerant 5/ states 
that participate in metallic bonding^. Alternatively, the 
heavy actinides (Am to Es) have larger atomic volumes, 
high-symmetry crystal structures, and relatively local- 
ized /-electrons. Applying pressure to the heavy ac- 
tinides results in a series of crystallographic phase tran- 
sitions, and the respective phases often have significantly 
different volumes'^''*. Transitions of this sort are often 
referred to as "volume collapse" transitions. Given that 
the application of ample pressure to any system of local- 
ized electrons will eventually cause a delocalization tran- 
sition, understanding what role the electronic delocaliza- 
tion transition may play in the volume collapse transition 
has been and continues to be an active area of study^'^. 

Plutonium is considered to be the dividing line of ac- 
tinide series, with the a- and 5- phases associated with 
light and heavy behavior, respectively. This dual nature 
of Pu, along with an enormous 25% volume collapse for 
the (5 — > a transition, has made Pu the most interesting 
element among the 5/ compounds for basic theoretical 
research over the past 50 years^'^. 

The actinides are among the most complicated classes 
of materials in terms of understanding electronic corre- 
lations given the presence of s, p, d, and / electrons 
near the Fermi surface and the unusual behavior ob- 
served in experiment. Broad discussion in the literature 
was devoted to the following topics: abrupt change in 
volume and bulk modulus^; unique crystal structures^; 
partial localization of /-electrons^*', Mott transition^^'^^; 
paramagnetism in light actinides and formation of mag- 
netic moments in heavier actinides (starting from Cm)®. 
For this purpose numerous Ab Initio electronic struc- 
ture calculations have been performed for the actinides. 
The techniques such as LBA^^'^^, GGA^^, SIC-LSDi^, 



LDA-hU""i9 and DMFT-based approachcs^^-^a j^^ve 
been implemented^'^. 



B. Actinide Hamiltonian 

The model Hamiltonian for the actinides can be writ- 
ten as: 

H = J2 ^^^-p + c-c.) + E iLi3fLf0j (1) 

where fai is the annihilation operator for 5/-electron in 
state a = \j,jz) at site i, and Cai is annihilation oper- 
ator for conduction electrons in the state a = \n,j',jg). 
This model can be understood as a periodic Anderson 
model in which additional direct hopping is allowed be- 
tween the correlated states. Equivalcntly, this model 
can be thought of as a Hubbard model with additional 
uncorrelated states that hybridize with the correlated 
states. Therefore, the model for the actinides contains 
the physics of both the periodic Anderson model and the 
Hubbard model. In the Hubbard model competes with 
U to determine the degree of localization of the electrons, 
while in the periodic Anderson model V competes with 
U. In the model of the actinides and V cooperatively 
compete with U, and the relative magnitudes of and V 
will determine the degree of Hubbard-like and Anderson- 
like contributions to the itineracy of the /-electrons. The 
main focus of this study is to determine the relative im- 
portance of and V across the actinide series. This is a 
first step towards a detailed understanding of the quanti- 
tative aspects of the localization-delocalization transition 
in the actinide series. 

Hamiltonians of the form described in Eq. (1) con- 
taining heavy and light electrons have appeared in vari- 
ous contexts in condensed matter physics. To deal with 
this complexity, this model is often reduced to a simpler 
model by eliminating the light electrons (i.e. spd states) 
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to obtain an effective Hubbard model having only the 
heavy (i.e. /) states. The reduced model only describes 
the bands within a narrow energy window around the 
Fermi energy (i.e. ^leV). The new renormalized / hop- 
pings have contributions from both the original / hop- 
ping in addition to the spd hopping. Additionally, many 
new interactions terms arc generated, but these are usu- 
ally ignored and only the on-site coulomb repulsion is 
retained as in the original hamiltonian. This procedure 
of one-electron " downfolding" has been used extensively 
to study electronic correlations in the transition-metal 
oxides, and there has been success in describing the pho- 
toemission spectra in this manner. In the context of the 
actinides, it is not clear that this approach is justified. 
Perhaps the most important issue is the nature of the 
localization transition in the actinidc series. In the Hub- 
bard model, when U is sufficiently large the effective / 
states will be localized and the system will be insulating. 
In the actinide model, when U is sufficiently large the / 
states will be localized, but the system will not neces- 
sarily be an insulator. It is possible that the spd states 
may still form a Fermi surface and give rise to a metal- 
lic state. Therefore, the actinide model and the effective 
Hubbard model differ even at a qualitative level in cer- 
tain regimes. Some aspects of localization-delocalization 
in the actinide model and the Hubbard model treated by 
DMFT are very similar at intermediate temperatures^'* 
(for example they both exhibit a line of first order phase 
transitions ending at a second order point), but there arc 
significant differences at very low temperatures when hy- 
bridization becomes a relevant perturbation suppressing 
Mott transition^^. Furthermore, the behavior at large U 
and high temperatures should be quite different in the 
two models, due to the presence of the broad metallic 
spd bands in the actinide model. Due to these consider- 
ations, we pursue the actinide model which should be a 
more accurate representations of the actinides given that 
the Hubbard model is an approximate reduction of the 
actinide model. 

In general, the parameters V and t-^ depend on the 
choice of basis set and therefore are not unique. This, 
of course, does not affect the band structure which is 
basis independent, but becomes an important in the 
context of approximate many body treatments (such as 
DMFT) which include only local Coulomb on the / or- 
bitals. For this purposes it is clearly advantageous to 
set up a Hamiltonian in an orthogonal basis where the 
/ electrons are highly localized. Hence, the secondary 
objective of this work is to determine a good basis set 
for setting up models of the actinides. For earlier tight- 
binding parametrization for actinides see Ref.^*'^^'^^. 



C. Motivation for this work 

The idea of a localization-delocalization transition in 
the actinides was brought forward by Johansson**. Jo- 
hansson based this idea on an empirical comparison of 



TABLE I: Choice of basis. 



Bare LMTO 

Lowdin transform 


Screened LMTO 

Lowdin transform 


Bare LMTO 

projective basis 


Screened LMTO 

projective basis 



the canonical 5/-bandwidth with the estimates of the 
Coulomb interaction in the form of a Hubbard U, and 
therefore he is starting with the assumption of a Hub- 
bard model to represent the actinides. As a result, the 
localization-delocalization transition is designated as a 
Mott transition in his paper, but it should be realized 
that this is a consequence of starting with an assumption 
of a Hubbard model. For a later elaboration of these ideas 
in the context of the a-S transition in Pu see Rcf.'^^'^^. 

The important role of d — / hybridization in actinide 
metals and alloys was stressed in the early work of Jullien 
et al.'^^''^* who considered models similar to Eq. (1). 

In this paper, we reconsider the issue of the descrip- 
tion of the localization-delocalization transition in the 
actinide series from a perspective which is motivated 
by recent DMFT and LDA-hDMFT works20-22^ ^hese 
works utilize the actinide hamiltonian (i.e. both / 
and spd states are included), and have provided fur- 
ther demonstration of the hypothesis that a localization- 
delocalization transition takes place across the actinide 
scries. However, the relative importance of and V (i.e. 
Hubbard- like vs. Anderson-like contributions) were not 
explicitly examined in these studies. 



II. ORBITALS AND BASIS 

A. Basis set dependence issue 

While the issue of representing the Kohn-Sham hamil- 
tonian in different basis sets has been a subject of numer- 
ous studies, the dependence of the results of correlated 
electronic structure methods such as LDA+DMFT on 
the choice of correlated orbitals is only beginning to be 
cxplored'^^ . 

In this study, we investigate the role of the choice of 
the correlated / orbital. We first take the / electron 
orbital as the / element of the LMTO basis, both in 
the bare and screened rcprcscntations'^'^''^^. The LMTO 
basis is non-orthogonal and therefore must be orthog- 
onalized in order to avoid the complications of solving 
the many-body problem in a non-orthogonal basis. As 
we will show in this study, the method of orthogonal- 
ization has a large influence on the results. We utilize 
both the Lowdin orthogonalization'^^ and the projective 
orthogonalization'^'^'''^ that was used in earlier implemen- 
tations of LDA-fDMFT. This effectively results in four 
different constructions of / orbitals, listed in Tabic I. 
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B. Bare and Screened LMTO within ASA scheme 

The basis set of linear mufSn-tin orbitals (LM- 
TOs) has been extensively used in electronic structure 
calculations"^^'"^^. Within the atomic sphere approxima- 
tion (ASA), LMTO is a minimal and efficient basis set 
with one basis function per site / and quantum pair 
L = {l,m). Although the LMTO method is physically 
transparent, the constructed basis is non-orthogonal. 

Below we sketch the derivation of the bare and screened 
LMTO basis set within the ASA. The construction of 
the bare LMTOs starts with so called envelope 

function^^, which is a decaying solution of the Laplace 
equation centered at the site /: 



KLiri) 



Ki{ri)YL{vi) = ( - 



i+i 



Ydh), (2) 



here r/ = r — R/, unit vector f/ indicates the direction 
of r/, Yl{ti) is a spherical function, and w is a scaling 
parameter associated with the linear size of unit cell. 

In any atomic sphere other than J, K]^{rj) can be rep- 
resented as: 



KL{ri) = -J2SiLj'L'JL'{ri'). 



(3) 



The function Jl{yi) = {ri /wYYLiii) stands for the reg- 
ular solutions of Laplace equation, and Sil.i'L' are struc- 
ture constants. 

Inside each atomic sphere we construct a linear com- 
bination of the solution ^/^(r/) of the Schrodinger equa- 
tion and its first derivative with respect to energy (/)j^(r/) 
at some fixed energy E^. 

The final step is to smoothly match the boundary con- 
ditions at the surface of sphere /: 



$f (r,) EE Af^cj^iLiri) + BfL<Pj^{vi) ^ Kl{vi) (4) 
and at the surface of sphere /' for all /' ^ I: 

(5) 

With the array of constants A and B determined from 
(4) - (5) we conclude the construction of bare LMTO 
basis function: 



X/L(r/) 



Xi(r/), r £ Interstitial. 



The Fourier transform of the LMTOs with respect to 
R/ — R/' gives: 



XkL(r) 



^fW-EL'*i'(r)^kLL', \r\<RMT, 
Zi^e^'^^KLir-R), \r\>RMT- 



(7) 



The standard LMTO method outlined above yields 
long-range orbitals. The concept of a screened LMTO 



was created to overcome the non-locality of the bare 
LMTO basis set^*. The method is based upon the idea 
of localizing the LMTOs by screening with multipoles 
added on the neighboring spheres. Namely, to each reg- 
ular solution of Laplace equation we add — Q/l of the 
irregular solution: 



J£iri)^JLiri)-aiLKL{ri). 



(8) 



The condition that on-site Laplace solution should not 
change leads to the Dyson-like equation for the screened 
structure constants: 



",a', 



(9) 



where matrix index a refers to the pair (/, L) and implies 
summation over repeated indices. The matrices aa = ai 
are diagonal for each I. In our calculations the choice 
of a's was as follows: as = 5.5166, ap = 0.5242, ad = 
0.1382 and a/ 0.0355. 

The screened and bare envelope functions are related 
by the transformation Ua'.a introduced in (9): 

Ktivj) = J2 KL'irr)[Si'L'.iL + a^, l' Sf, ^, j^], (10) 

I'L' 

Or in matrix notations: 

where Ka = KL{ri) and = Kl{ri). 

With the definitions (8), (9) and (10) the construction 
of screened LMTOs proceeds exactly in the same way as 
in the case of the bare LMTOs. Namely, we construct 
new linear combinations 

inside the sphere / by matching smoothly X£(r/) on its 
surface. Also, we construct new linear combination 

inside sphere /' by matching smoothly J£, (r//) on its 
surface for each /' ^ I. Thus, we arrive to the definition 
of the screened LMTO: 



(6) 



X?L(r/) 



$f"(r,), rjeSi, 

-Ei'L' Shj'L'^inri'), rj, G Sril ^ /'),(11) 

K1{ri), r G Interstitial. 



The Fourier transform of the screened LMTOs with 
respect to R/ — R// gives: 



*f " W - El' <i>i?(r) |r| < RmT: 

j:^e'^^K2{r-K), \r\ > Rmt- 



(12) 
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The Hamiltonian and overlap matrices in screened and 
bare LMTO representations {O , H, and O", if" respec- 
tively) are related through the transformation U intro- 
duced in (9): 



H°' = U^'HU 
O" = U'^OU. 



(13) 
(14) 



Having constructed the basis, one has to solve the gen- 
eralized eigenvalue problem: 



(if(k)-e,(k)O(k))^',(k)==0. 



(15) 



As described above, it is necessary to transform to an or- 
thogonal basis when performing many-body calculations, 
such as DMFT, in order to avoid the difficulties associ- 
ated with a non-orthogonal basis. 



C. Lowdin orthogonalization 

Lowdin orthogonalization^^ is a straightforward or- 
thogonalization of the Hamiltonian which uses no infor- 
mation from the basis set: 



H(k) 



1 



W(k)- 



1 



(16) 



As will be shown below, this orthogonalization proce- 
dure may lead to a further mixing of L characters among 
the LMTOs and hence unphysical results. 



D. Projective orthogonalization 



In each L-subspace the overlapping matrix is: 

Ok = o^""^ - 5^o(^^) - o(^'^)5k + Slo^'-'^S^. (19) 

In order to find the transformation to the orthonormal 
base we must represent 0(k) as the square of a matrix. 
As we show below, the most intelligent choice would be: 



0(k) « (7^- J5k)t(H- J5k) 



(20) 



for each i-subspace. Here Ti, and J are diagonal matrices 
proportional to unity in each subspace of definite L just 
like the overlaps o'^^^^ defined above. 

The above equation can not be made exact because 
the overlap numbers are obtained by integration over the 
radial part of wave functions. However, in most cases the 
overlap numbers can become very close to their approxi- 
mations: 



{HH) 
(JH) 

JJJ) 



ril riU 



(21) 



For each L we have three independent equations for 
two unknowns. An approximate solution can be found 
by minimizing the following function: 



, n^Jil' + |oJ""' - Ji*Ji\^ = mm. 

The desired transformation to the new base is: 



(HJ) 



ntHii^ + \o\-'"' - ji*ni\ 
2 I iJ-'-n 



(22) 



A physically motivated orthogonalization procedure is 
to find a basis where each function contains the maxi- 
mum amount of a particular L character. This approach 
proposed by K. Haule was used in earlier LDA-I-DMFT 
studies of Cerium and Plutonium^^. This basis has an 
important advantage, being that the "/ electron" in this 
basis has maximal / character. Mathematically, the non- 
interacting spectral function of the / electron Green's 
function in this basis agrees with the LDA projected 
density-of-states having / character, as shown in Ap- 
pendix A. This allowed us to identify the / occupation 
in this basis set with the occupation numbers inferred 
from EELS and X-Ray absorption which are sensitive to 
angular momentum selection rules^^. 

Here we follow Ref.'^^. It is straightforward using (7) 
and (14) to show that overlapping matrix within MT- 
sphere can be represented as^^: 



(JH) 



o|f'^^kL,L.(17) 



(J. J) n 
0„ OkL'Lo 



The quantities of^, of^ , of"^ and of'^ are numbers in 
each /-subspace. For A and B representing H or J: 



(18) 



Finally: 



Tk = (H - JSkY 



(23) 



(24) 



III. RESULTS 

We perform relativistic, spin-restricted LDA calcula- 
tions within the ASA scheme. 7s, 6p, 6d and 5/-orbitals 
were chosen to represent valence states, and lO'^k-points 
were used in the first Brillouin zone. The same type 
of calculations were carried out for 4 different materials, 
picked to evenly represent actinidc series: U, a-Pu, S- 
Pu and Cm II {fee phase of curium). For simplicity, we 
used the fee crystal structure for each element. The lat- 
tice parameters listed in Table II were chosen to match 
the experimentally measured volumes for corresponding 
phases in case of Pu and Cm II. For U we use the equi- 
librium volume predicted within GGA calculations. 

In Figure 1 we present LDA band structure of Cm 
II with the projections of the 5/- and 6d- characters. 
The overwhelming contribution of / character within a 
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TABLE II: Lattice parameters(in rA). 



a-U 


4.3378 


a-Pu 


4.3074 


S-Pu 


4.6400 


Cm II 


4.9726 




leV window around Fermi level suggests the conclusion 
that the low-energy physics of actinides is completely 
controlled by / — / bonding. As we show below this 
intuitive interpretation turns out to be misleading and 
Hubbard model alone can not be considered as the low- 
energy Hamiltonian for actinides. One has to account for 
presence of spd-characters at the Fermi level through the 
hybridization. Moreover, the hybridization energy scale 
in actinides turns out to be larger than the average / — / 
hopping. 



W X K 



FIG. 1: (Color online) Band structure of Cm with indicated 
contribution of 5/ (light shade, red online) and 6d (dark 
shade, green online) characters. The layout is chosen to show 
c?-contributions over /-contribution. Inset: Brillouin zone of 
fee structure with indicated high symmetry directions. 



A. Determining a robust basis for the actinides 



In order to determine the optimum basis, we need to 
define a criteria to judge the different bases. When per- 
forming DMFT calculations, one accounts only for a sub- 
set of local electronic correlations (those on the / orbital) . 
Therefore, from the perspective of DMFT it is best to 
have / orbitals with the largest on-site Coulomb repul- 
sion A simpler criteria, in the same spirit, is to 
search for the smallest value of . We first investigate 

in Cm for the four different basis sets in Table I. The 
hybridization V in the Hamiltonian (1) may be set to 
zero. What remains are the two blocks Hf and Hspd 
which are now completely decoupled. The Hamiltonian 
may now be diagonalized resulting in distinct spd and 
/ bands, and any dispersion of the / bands is due to 

. We begin by analyzing the bare LMTOs orthogo- 
nalized with the Lowdin procedure (see top left panel 
Figure 2). Some / bands have a dispersion greater than 
1.5eV which is unfavorable. Using the bare LMTOs or- 
thogonalized with the projective procedure, the / bands 
are far more narrow with a width of less than OAeV (see 
left bottom panel of Figure 2). In this case the two sets 
of bands can be identified as S" = | and S* = |. The 
Lowdin orthogonalization mixes the spd states into the 
/ states which causes a larger dispersion and a mixing 
of / bands between the S = ^ and 5=1 states. Alter- 
natively, the projective orthogonalization minimizes the 



amount of spd character in the / states which results in 
weakly dispersing / states. 

The same exercise can be performed using the screened 
LMTOs (see right top and bottom panels of Figure 2). 
In this case, both the Lowdin and the projective orthogo- 
nalization produce nearly identical results to the projec- 
tive orthogonalization of the bare LMTOs. The screened 
LMTOs are insensitive to the method of orthogonaliza- 
tion due to the fact that orbitals are already well localized 
with a well-defined character. In conclusion, one may use 
bare LMTOs orthogonalized with the projective proce- 
dure or screened LMTOs orthogonalized in an arbitrary 
manner as a robust basis for the actinides. 
B. Decomposition of the actinide band structures 

Having established a sensible basis for the actinides, 
we choose to proceed with projective orthogonalization of 
bare LMTOs. It is instructive to zero the hybridization 
V of the Hamiltonian for U, a-Pu, S-Pu, and Cm, and to 
compare the full band structure with the spd and / bands 
(see Figures 3 and 4). The same generic behavior can be 
seen in all four systems. The spd bands have a strong 
dispersion and cross the Fermi energy in all cases, and 
the / bands are relatively narrow. The fact that the spd 
bands cross the Fermi energy in all cases is a critical point 
which indicates that there will be spd states at the Fermi 
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FIG. 2: (Color online) Basis difference for fee Curium. In all panels dashed grey lines represent the LDA band structure , solid 
light shade lines (green online) represent bands for the block Hs^d , and solid dark shade line (red online) represent bands for 
the block Hf. 



energy even if the / states become completely localized. 
When the hybridization Y is switched on, the / and spd 
bands interact via V and mix. Therefore the strength of 
V can qualitatively be seen as the difference between the 
full DFT bands and the f+spd bands. The / bands are 
relatively wide for Uranium and become increasingly nar- 
row as Curium is approached. The spd bands follow the 
same general trend, but the relative changes are smaller. 
The values of V and will be quantified below. 



C. Quantitative analysis of V and 

In order to quantify V and for the different actinides, 
we introduce an average V and so each actinidc may 
be characterized by two numbers. 

First, we recall that the Hamiltonian (1) consist of four 
blocks: 



Hik) 



[ H/(k) 



(25) 



Then the average strength of the hybridization per 
band is defined as follows: 



V 



(26) 



where H{k) stands for hamiltonian (25) with Hspd{k) = 
Hf{k) ~ 0, Nf ~ 14 stands for number of /-bands, and 
(•■•) = tI; Ek ■•• • The definition (26) was chosen to 
match hybridization V of standard Anderson model in 
two- band limit. 

The average value of is defined as follows: 



tf = ^[Tr{{Hf{kr) ^ {Hf{k)r)] 



1/2 



(27) 



and matches of the canonical Hubbard model in the 
limit of one-band model. 

Table III lists the calculated values of the average 
hybridization V and and the average energy for the 
j = 5/2 and j = 7/2 levels of the / manifold relative to 
the Fermi energy. The averages are generally the same 
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FIG. 3: (Color online) Band structure of a-U (top) and Cm 
II (bottom). Grey dashed lines represent LDA bands, solid 
light shade lines (green online) represent bands of Hspd, and 
solid dark shade lines (red online) represent bands of Hj. 



FIG. 4: (Color online) Band structure of q-Pu (top) and 5- 
Pu(bottom). Grey dashed lines represent LDA bands, solid 
light shade lines (green online) represent bands of Hspd, and 
solid dark shade lines (red online) represent bands of 



TABLE III: Quantitative characteristics for actinide series (in 
eV). 





V 


tf 


V/tl 


£5/2 - M 


£7/2 - ^J■ 


Bare LMTO 


Q-U 


0.483 


0.188 


2.569 


0.442 


1.353 


a-Pu 


0.423 


0.146 


2.897 


-0.180 


0.971 


5-Pu 


0.305 


0.099 


3.081 


-0.129 


1.008 


Cm II 


0.189 


0.050 


3.780 


-1.152 


0.238 


Screened LMTO 


a-U 


0.490 


0.188 


2.606 


0.444 


1.355 


a-Pu 


0.429 


0.146 


2.938 


-0.178 


0.973 


5-Pu 


0.309 


0.098 


3.153 


-0.128 


1.009 


Cm II 


0.192 


0.050 


3.840 


-1.151 


0.238 



for the bare and screened LMTOs, with the exception of 
the average hybridization being slightly larger in the case 
of screened LMTOs. 



These results are displayed graphically in Figure 5. In 



TABLE IV: Nearest Neighbors contributions to V and tf (in 
eV). 





v„ 


tf 

''71 




Q-U 


0.371 


0.172 


2.157 


Q-Pu 


0.324 


0.134 


2.418 


<S-Pu 


0.232 


0.090 


2.578 


Cm II 


0.143 


0.045 


3.178 



all cases, V is significantly greater than . As one moves 
along actinides series from U to Cm tf decreases as much 
as four times. The average value of hybridization V also 
decreases but at a slower rate, as indicated by the in- 
set plot of the ratio of V and t^ . The strong decrease 
in V and tf will both contribute to the localization of 
the / states. In order to determine if the localization 
could be predominantly assigned to either Mott or An- 
derson character, explicit many-body calculations such 
as DMFT would need to be performed. 

In order to provide further insight into the degree of 
locality of the basis, it is instructive to determine the frac- 
tion of Vn and tn which arise solely from nearest-neighbor 
hopping. The corresponding values are presented in Ta- 




a-U a-Pu 6-Pu Cm II 



FIG. 5: (Color online) Histogram represents average hy- 
bridization (first bar for each element, blue online) and aver- 
age f — f hopping (second bar for each element, red online) 
as functions of atomic number. The shadow bars show V and 
ff for the original Hamiltonian and bright bars represent val- 
ues Vn and til arising from nearest neighbors contributions 
only. In inset: the ratio V /tf (squares) as function of atomic 
number. The ratio of nearest neighbor contributions is repre- 
sented by circles. 

blc IV and can also be seen in Figure 5. First nearest 
neighbors contribute ?» 75% to V and w 90% to the . 

The ratio of is also given for the nearest-neighbor 

contribution, and the shape and slopes of the two respec- 
tive curves are very similar. This analysis indicates that 
nearest-neighbor hopping in real space accounts for most 
of the relevant one-electron physics. 



D. Real space analysis of band structure 

In the above analysis it was shown that nearest- 
neighbor hopping accounts for a strong majority of V 
and tJ . Therefore, it is suggestive that the one-electron 
bands can be reproduced with relatively short ranged 
hoppings and V. In order to determine the degree of 
locality, we plot the band structure as a function of the 
number of neighbors for the t-l^ and V hoppings (see Fig- 
ures 6 and 7). Results are given for zero, one, and fourth 
neighbor hopping. The hopping f'P'^ are not truncated as 
it is clear they will definitely have relatively long-range 
hoppings. The generic results are similar for all four ma- 
terials. The case with zero neighbors results in flat bands 
having the on-site energy of each orbital. When first 
nearest-neighbors are included the resulting bands are an 
excellent approximation to the full band structure. In- 
cluding up to four nearest-neighbors yields nearly perfect 
agreement. Cm has better agreement than U for a given 
number of nearest-neighbors, and this reflects the larger 
degree of localization in the late actinides as compared to 
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5-Plutonium 
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FIG. 6: (Color online) Band structure of Curium (top) and 5- 
Pu (bottom) when N nearest neighbors are taken into account 
for /-orbitals (red solid line) is compared to original LDA 
bands (dashed grey line). The band structures are plotted for 
T-X direction. 



the early actinides. In conclusion, the band structure of 
the actinides is dominated by nearest-neighbor hopping 
when using an appropriate basis. 

E. Conclusion 

In summary, a one-electron analysis of band structure 
of the actinides was presented. We demonstrated that 
bare LMTOs orthogonalized with the projective method 
and screened LMTOs arc robust bases, in the sense that 
they give rise to / orbitals with minimal hopping. Analy- 
sis of the Hamiltonian in these bases yielded a number of 
interesting results. When switching off the hybridization 
V, it was shown that the spd states cross the Fermi en- 
ergy and hence will be present at the Fermi energy even 
if the / electrons become localized. Our description is in 
reasonable agreement with the earlier work of Harrison^^. 
In particular, the matrix elements of spin-orbit coupling 
indeed have atomic-like nature and the hybridization is 
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FIG. 7: (Color online) Band structure of q-Pu (top) and a-U 
(bottom) when N nearest neighbors are taken into account for 
/-orbitals (red solid line) is compared to original LDA bands 
(dashed grey line). The band structures are plotted for V-X 
direction. 



much larger than the direct / — / hopping^^ . However the 
spd bands are not simple plane waves and the hybridiza- 
tion matrix clement does not have a simple k-dependence 
proportional to an ^ = 2 spherical harmonic. 

Evaluation of the average hybridization V and average 
/ — / hopping as a function of the actinides showed 
that both quantities decrease strongly. The quantity 
decreased faster than V ^ but V was larger in all ac- 
tinides. Hence, the Anderson model of the localization- 
delocalization transition, rather than a multiorbital Hub- 



^w(^) -^-Y. {(^L - Gk.)o^^ - [5k(GL - Gk.)] 

k 



where the momentum dependent Green's function is 

Gk^-(Ok(^ + Ai)-i?k-Sk^)"' (A4) 



bard model is needed to describe the physics of the ac- 
tinides once explicit many-body calculations are added. 
This is the point of view taken in recent DMFT work^", 
and no further reduction to a model containing only f 
bands seems possible. Finally, a real-space analysis of 
the band structure demonstrated that nearest-neighbor 
hopping accounts for most of the band structure in the 
basis used in this study, thus providing a tight binding 
fitting of the bands of the actinides that can be useful in 
further studies. 
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APPENDIX A: GREEN'S FUNCTION IN THE 
PROJECTIVE BASE 

In the DMFT approach one needs to choose a set 
of localized Wannier states in which correlations are 
strongest. In the context of actinides, the orbitals with 
the largest component of 5/ character are the appropri- 
ate set of orbitals. 

For many-body calculations, it is convenient if the set 
of localized orbitals is orthogonal. In this case, it is de- 
sired that the local Green's function is connected to the 
partial density of states by the usual relation 

-DLm'(^) ~ TT-XGIoc - Gloc)lni,lm'- (Al) 
ZTTl 

In another words, the localized set of orbitals need to 
give rise to the 5/ spectra defined by 

^wM = / ^|^l^*™(fi)(Gt(ri,r2)-G(ri,r2))y,™,(f2) 

with I = 3 for actinides. Only in this case, the number of 
/ electrons (or the valence of the material) is connected 
to the impurity / count, as obtained in the DMFT cal- 
culation. 

Using the LMTO basis set Eq. (7), the partial density 
of states becomes 



UGL G^^)Sly" + [5k(GL - G^^)Sl]o'''} (A3) 
I 

and overlap numbers o'^^ are defined by Eq. (18). 
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The projective orthogonalization Eq. (24) leads to the = - JS\^)G\^^{H - JS'k)'''. (A6) 

foUowing Green's function 

^ The local spectral function in this new base therefore 

Gkc. = {Tl[O^,{L0 + /i) - Fk - Sk^]Tk)-^ (A5) becomes 

I 



^^j^l^Zoc ~ Gloc)lra,hn' - ^ { I^'L ^ G\^u\lL''Hi'Hi - [S\^{G\^^ - G]^^)\lL''Hi Jl 

k 

-[(gL - G^^)Si]Lvm + [5k(GL - G^^)Si]LvJ;Ji] (A7) 

I 

which is equivalent to the partial density of state Eq. (A3) accuracy (better than 1%) therefore the relation between 

provided the condition Eq. (21) is satisfied. Extensive the partial density of states and local Green's function 

experience shows that in the case of localized d and / Eq. (Al) is also satisfied to high accuracy, 
orbitals, the condition is always satisfied to very high 
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